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Abstract 

The branching fractions of the decay B + — > ppK + for different intermediate states 
are measured using data, corresponding to an integrated luminosity of l.Ofb -1 , 
collected by the LHCb experiment. The total branching fraction, its charmless 
component (M p p < 2.85GeV/c 2 ) and the branching fractions via the resonant cc 
states rj c (lS) and ip(2S) relative to the decay via a J/ifj intermediate state are 

B(B+ -> P pK+) total _ 4 gi ± Q ig (gtat) ± Q M (gygt)i 

: 2.02 ± 0.10 (stat) ± 0.08 (syst), 
: 0.578 ± 0.035 (stat) ± 0.027 (syst), 
: 0.080 ± 0.012 (stat) ± 0.009 (syst). 
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Upper limits on the i? + branching fractions into the rj c (2S) meson and into the 
charmonium-like states X(3872) and X(3915) are also obtained. 
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1 Introduction 



The B + — > ppK + decajj^] offers a clean environment to study cc states and char monium- like 
mesons that decay to pp and excited A baryons that decay to pK + , and to search for 
glueballs or exotic states. The presence of pp in the final state allows intermediate states 
of any quantum numbers to be studied and the existence of the charged kaon in the final 
state significantly enhances the signal to background ratio in the selection procedure. 
Measurements of intermediate charmonium-like states, such as the X(3872), are important 
to clarify their nature jlj[2] and to determine their partial width to pp, which is crucial to 
predict the production rate of these states in dedicated experiments (31. BaBar and Belle 
have previously measured the B + — > ppK + branching fraction, including contributions 
from the J/ip and rj c {lS) intermediate states |4j[5j . The data sample, corresponding to 
an integrated luminosity of l.Ofb -1 , collected by LHCb at y/s = 7TeV allows the study 
of substructures in the B + — > ppK + decays with a sample ten times larger than those 
available at previous experiments. 

In this paper we report measurements of the ratios of branching fractions 

B{B+ -» mode -» ppK+) 
n{m ° de) = B(B+^J^K+^ppK+y (1) 

where "mode" corresponds to the intermediate T] C (1S), i[)(2S), T] c (2S), XdoO-P)i h c (lP), 
X(3872) or X(3915) states, together with a kaon. 



2 Detector and software 

The LHCb detector |6| is a single-arm forward spectrometer covering the pseudorapidity 
range 2 < rj < 5, designed for the study of particles containing b or c quarks. The detector 
includes a high precision tracking system consisting of a silicon-strip vertex detector 
surrounding the pp interaction region, a large-area silicon-strip detector located upstream 
of a dipole magnet with a bending power of about 4 Tm, and three stations of silicon-strip 
detectors and straw drift-tubes placed downstream. The combined tracking system has 
momentum (p) resolution Ap/p that varies from 0.4% at 5GeV/c to 0.6% at lOOGeV/c, 
and impact parameter resolution of 20 |xm for tracks with high transverse momentum 
(pr)- Charged hadrons are identified using two ring- imaging Cherenkov (RICH) detectors. 
Photon, electron and hadron candidates are identified by a calorimeter system consisting of 
scintillating-pad and pre-shower detectors, an electromagnetic calorimeter and a hadronic 
calorimeter. Muons are identified by a system composed of alternating layers of iron and 
multiwire proportional chambers. 

The trigger 17) consists of a hardware stage, based on information from the calorimeter 
and muon systems, followed by a software stage where candidates are fully reconstructed. 
The hardware trigger selects hadrons with high transverse energy in the calorimeter. The 
software trigger requires a two-, three- or four-track secondary vertex with a high p^ sum 

1 The inclusion of charge-conjugate modes is implied throughout the paper. 
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of the tracks and a significant displacement from the primary pp interaction vertices (PVs). 
At least one track should have px > 1.7GeV/c and impact parameter (IP) y 2 with respect 
to the primary interaction greater than 16. The IP y 2 is defined as the difference between 
the x 2 °f the PV reconstructed with and without the considered track. A multivariate 
algorithm is used for the identification of secondary vertices consistent with the decay of a 
b hadron. 

Simulated B + — > ppK + decays, generated uniformly in phase space, are used to 
optimize the signal selection and to evaluate the ratio of the efficiencies for each considered 
channel with respect to the J/ip channel. Separate samples of B + — > J/ipK + — > ppK + 
and B + — > r] c {lS)K + — > ppK + decays, generated with the known angular distributions, 
are used to check the dependence of the efficiency ratio on the angular distribution. In 
the simulation, pp collisions are generated using Pythia 6.4 |8| with a specific LHCb 



configuration 91. Decays of hadronic particles are described by EvtGen 10 in which final 



state radiation is generated by Photos 11 . The interaction of the generated particles 
with the detector and its response are implemented using the Geant4 toolkit 12 as 



described in Ref. 13 1. 



3 Candidate selection 

Candidate B + — > ppK + decays are reconstructed from any combination of three charged 
tracks with total charge of +1. The final state particles are required to have a track fit 
with a x 2 /ndf < 3 where ndf is the number of degrees of freedom. They must also have 
p > 1500MeV/c, pt > lOOMeV/c, and IP y 2 > 1 with respect to any primary vertex 
in the event. Particle identification (PID) requirements, based on the RICH detector 
information, are applied to p and p candidates. The discriminating variables between 
different particle hypotheses (tt, K, p) are the differences between log-likelihood values 
Aln£ Q(9 under particle hypotheses a and (3, respectively. The p and p candidates are 
required to have A\nC p7T > —5. The reconstructed B + candidates are required to have an 
invariant mass in the range 5079 — 5579 MeV/c 2 . The asymmetric invariant mass range 
around the nominal B + mass is designed to select also B + — > ppir + candidates without any 
requirement on the PID of the kaon. The PV associated to each B + candidate is defined 
to be the one with the smallest IP y 2 . The B + candidate is required to have a vertex fit 
with a y 2 /ndf < 12 and a distance greater than 3 mm, a y 2 for the flight distance greater 
than 500, and an IP y 2 < 10 with respect to the associated PV. The maximum distance of 
closest approach between daughter tracks has to be less than 0.2 mm. The angle between 
the reconstructed momentum of the B + candidate and the B + flight direction (Or) is 
required to have cos^a > 0.99998. 

The reconstructed candidates that meet the above criteria are filtered using a boosted 



decision tree (BDT) algorithm 14 . The BDT is trained with a sample of simulated 
B + — > ppK + signal candidates and a background sample of data candidates taken from 
the invariant mass sidebands in the ranges 5080 — 5220MeV/c 2 and 5340 — 5480 MeV/c 2 . 
The variables used by the BDT to discriminate between signal and background candidates 
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BDT response 



Figure 1: Distribution of the BDT algorithm response evaluated for background candidates 
from the data sidebands (red), and signal candidates from simulation (blue). The black dotted 
line indicates the chosen BDT response value. 

are: the of each reconstructed track; the sum of the daughters' p?', the sum of the IP \ 2 
of the three daughter tracks with respect to the primary vertex; the IP of the daughter, 
with the highest pt ; with respect to the primary vertex; the number of daughters with 
Pt > 900GeV/c; the maximum distance of closest approach between any two of the B + 
daughter particles; the IP of the B + candidate with respect to the primary vertex; the 
distance between primary and secondary vertices; the 6$ angle; the xV n df of the secondary 

vertex; a pointing variable defined as D . p a S1 ^ , where P is the total momentum of the 

three-particle final state, 9 is the angle between the direction of the sum of the daughter's 
momentum and the direction of the flight distance of the B + and ^ pr,i is the sum of the 
transverse momenta of the daughters; and the log likelihood difference for each daughter 
between the assumed PID hypothesis and the pion hypothesis. The selection criterion on 
the BDT response (Fig. [I]) is chosen in order to have a signal to background ratio of the 
order of unity. This corresponds to a BDT response value of —0.11. The efficiency of the 
BDT selection is greater than 92% with a background rejection greater than 86%. 

4 Signal yield determination 

The signal yield is determined from an unbinned extended maximum likelihood fit to 
the invariant mass of selected B + — > ppK + candidates, shown in Fig. [2^,). The signal 
component is parametrized as the sum of two Gaussian functions with the same mean 
and different widths. The background component is parametrized as a linear function. 
The signal yield of the charmless component is determined by performing the same fit 
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Figure 2: Invariant mass distribution of a) all selected B + — > ppK + candidates and b) candidates 
having M p p < 2.85GeV/c 2 . The points with error bars are the data and the solid lines are the 
result of the fit. The dotted lines represent the two Gaussian functions (red) and the dashed line 
the linear function (green) used to parametrize the signal and the background, respectively. The 
vertical lines indicate the signal region. The two plots below the mass distributions show the 
pulls. 



described above to the sample of B + — > ppK + candidates with M pp < 2.85GeV/c 2 , shown 
in Fig. [2)3). The B + mass and widths, evaluated with the invariant mass fits to all of 
the B + — > ppK + candidates, are compatible with the values obtained for the charmless 
component. 

The signal yields for the charmonium contributions, B + — > (cc)K + — > ppK + , are 
determined by fitting the pp invariant mass distribution of B + — > ppK + candidates within 
the B + mass signal window, \M ppK + — M B +\ < 50MeV/c 2 . Simulations show that no 
narrow structures are induced in the pp spectrum as kinematic reflections of possible 
B + — > pA — > ppK + intermediate states. 

An unbinned extended maximum likelihood fit to the pp invariant mass distribution, 
shown in Fig. [3J is performed over the mass range 2400 — 4500 MeV/c 2 . The signal 
components of the narrow resonances J/ip, ip(2S), h c (lP), and X(3872), whose natural 
widths are much smaller than the pp invariant mass resolution, are parametrized by 
Gaussian functions. The signal components for the ?7 C (15'), Xco(l-P); Vci^S), and X(3915) 
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Figure 3: Invariant mass distribution of the pp system for B + — > ppK + candidates within 
the B + mass signal window, \M(ppK + ) — M B + \ < 50MeV/c 2 . The dotted lines represent the 
Gaussian and Voigtian functions (red) and the dashed line the smooth function (green) used to 
parametrize the signal and the background, respectively. The bottom plot shows the pulls. 



are parametrized by Voigtian functions^] Since the pp invariant mass resolution is ap- 
proximately constant in the explored range, the resolution parameters for all resonances, 
except the i[)(2S), are fixed to the J/ip value (<Jj/^ = 8.9 ± 0.2 MeV/c 2 ). The background 
shape is parametrized as f{M) = e c ^ M+c ^ M2 where c\ and C2 are fit parameters. The J/ip 
and i^(2S) resolution parameters, the mass values of the rj c (lS), J/ip, and ip(2S) states, 
and the r) c (lS) natural width are left free in the fit. The masses and widths for the other 
signal components are fixed to the corresponding world averages [15] . The pp invariant 
mass resolution, determined by the fit to the ip(2S) is a^(2S) — 7.9 ± 1.7MeV/c 2 . 

The fit result is shown in Fig. |3j Figures [4] and [5] show the details of the fit result in 
the regions around the rj c (lS) and J/ip, t] c {2S) and i[)(2S), Xco(l-P) an d h c , and X(3872) 
and X(3915) resonances. Any bias introduced by the inaccurate description of the tails of 
the ?7 C (15'), J/ip and ip{2S) resonances is taken into account in the systematic uncertainty 
evaluation. 



2 A Voigtian function is the convolution of a Breit-Wigner function with a Gaussian distribution. 
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Figure 4: Invariant mass distribution of the pp system in the regions around a) the rj c (lS) 
and J/ifi and b) the r] c (2S) and ip(2S) states. The dotted lines represent the Gaussian and the 
Voigtian functions (red) and the dashed line the smooth function (green) used to parametrize 
the signal and the background, respectively. The two plots below the mass distribution show the 
pulls. 
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Table 1: Signal yields for the different channels and corresponding 95% CL upper limits for 
modes with less than 3a statistical significance. For the J/tp mode, the non-signal yield is 
subtracted. Uncertainties are statistical only. 



B + decay mode 


Signal yield 


Upper limit (95% CL) 


ppK+ [total] 


6951 ± 176 




ppK+ [M pP < 2.85 GeV/c 2 ] 


3238 ± 122 




J/ipK+ 


1458 ± 42 




Vc(lS)K+ 


856 ±46 




i){2S)K + 


107 ±16 




Vc (2S)K+ 


39 ±15 


< 65.4 


Xco(lP)K+ 


15 ±13 


< 38.1 


h c (lP)K+ 


21 ±11 


< 40.2 


X(3872)K+ 


-9±8 


< 10.3 


X(3915)K + 


13 ±17 


< 42.1 



The contribution of cc — > pp from processes other than B + — > ppK + decays, denoted 
as "non-signal" , is estimated from a fit to the B + mass in the sidebands 5130 — 5180 and 
5380 — 5430MeV/c 2 . Except for the J/ip mode, no evidence of a non-signal contribution 
is found. The non-signal contribution to the J/ip signal yield in the B + mass window is 
43 ± 11 candidates and is subtracted from the number of J /ip signal candidates. 

The signal yields, corrected for the non-signal contribution, are reported in Table [I] 
For the intermediate charmonium states r] c (2S), x c0 (lP), h c (lP), X(3872) and X(3915), 
there is no evidence of signal. The 95% CL upper limits on the number of candidates 
are shown in Table [T] and are determined from the likelihood profile integrating over the 
nuisance parameters. Since for the X(3872) the fitted signal yield is negative, the upper 
limit has been calculated integrating the likelihood only in the physical region of a signal 
yield greater than zero. 



5 Efficiency determination 

The ratio of branching fractions is calculated using 

, x B(B + -)• mode -)• ppK + ) N mode e^, 
K{mode) - B{B+ - Jf ^ K+ - p _ K+) - — x — (2) 

where A^ mo d e and Nj/^ are the signal yields for the given mode and the reference mode, 
B + — > J/tpK + — > ppK + , and e mo d e /ej/^ is the corresponding ratio of efficiencies. The 
efficiency is the product of the reconstruction, trigger, and selection efficiencies, and is 
estimated using simulated data samples. 

Since the track multiplicity distribution for simulated events differs from that observed 
in data, simulated candidates are assigned a weight so that the weighted distribution 
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Figure 6: Efficiency as a function of M p p for B + — > ppK + decays. The solid line represents the 
linear fit to the efficiency distribution; the dashed line is the point-by-point interpolation used to 
estimate the systematic uncertainty. 

reproduces the observed multiplicity distribution. The distributions of Aln£ Kn and Aln£ p7r 
for kaons and protons in data are obtained in bins of momentum, pseudorapidity and 
number of tracks from control samples of D* + — > £)°(— > K~tt + )tt + decays for kaons and 
A — > pn~ decays for protons, which are then used on a track-by-track basis to correct the 
simulation. The efficiency as a function of M p p is shown in Fig. [6j A linear fit to the 
efficiency distribution is performed and the efficiency ratios are determined based on the 
fit result. 

6 Systematic uncertainties 

The measurements of the relative branching fractions depend on the ratios of signal yields 
and efficiencies with respect to the reference mode. Since the final state is the same in 
all cases, most of the systematic uncertainties cancel. The systematic uncertainty on the 
efficiency ratio, in each region of pp invariant mass, is determined from the difference 
between the efficiency ratios calculated using the solid fitted line and the dashed point-by- 
point interpolation shown in Fig. [6j The uncertainty associated with the evaluation of 
the B + signal yield has been determined by varying the fit range by ±30MeV/c 2 , using a 
single Gaussian instead of a double Gaussian function to model the signal PDF, and using 
an exponential function to model the background. For each charmonium resonance the 
systematic uncertainty on the signal yield has been investigated by varying the B mass 
signal window by ±10MeV/c 2 , the signal and background shape parametrization and 
the subtraction of the cc contribution from the continuum. The systematic uncertainty 
associated with the parametrization of the signal tails of the J/if>, r] c (lS) and ip^S) 
resonances is taken into account by taking the difference between the number of candidates 




9 



Table 2: Relative systematic uncertainties (in %) on the relative branching fractions from 
different sources. The total systematic uncertainty is determined by adding the individual 
contributions in quadrature. 



Source 


ft(total) 


Tl{Mpp < 2.85 GeV/c 2 ) 


K( Vc (lS)) 


W(2S)) 


Efficiency ratio 


0.21 


0.5 


3.3 


4.8 


B + mass fit range 


0.16 


0.5 






Sig. and Bkg. shape 


2.5 


3.6 


1.8 


6.5 


B + mass window 


0.6 


0.6 


0.9 


3.8 


Non-signal component 






0.4 


5.1 


Signal tail param. 


1.0 


1.0 


1.2 


4.3 


Total 


2.8 


3.8 


4.1 


11.3 



Source 


n(Vc(2S)) 


ft(Xco(lP)) 


n(h c (ip)) 


ft(X(3872)) 


ft(X(3915)) 


Efficiency ratio 


4.4 


2.5 


3.4 


6.5 


7.0 


B + mass fit range 












Sig. and Bkg. shape 


3.9 


3.3 


14.3 


5.6 


10.1 


B + mass window 


11.3 


23.6 


23.6 


17.5 


7.5 


Non-signal component 












Signal tail param. 


1.0 


1.0 


1.0 


1.0 


1.0 


Total 


12.8 


24.0 


27.8 


19.5 


15.5 



in the observed distribution and the number of candidates calculated from the integral 
of the fit function in the range —6a to —2.5a. The systematic uncertainty associated 
with the selection procedure is estimated by changing the value of the BDT selection 
to —0.03, which retains 85% of the signal with a 30% background, and is found to be 
negligible. The contributions to the systematic uncertainties from the different sources are 
listed in Table |2j The total systematic uncertainty is determined by adding the individual 
contributions in quadrature. 



7 Results 

The results are summarized in Table [3] and the values of the product of branching fractions 
derived from our measurement using the world average values B(B + — > J/ipK + ) = 
(1.013 ± 0.034) x 10" 3 and B(J/i/j ->■ pp) = (2.17 ± 0.07) x 10~ 3 |15) are listed in Table 



The branching fractions obtained are compatible with the world average values 15 . The 
upper limit on B(B + — > Xco(^-P)K + — > ppK + ) is compatible with the world average 
B(B+ X co(lP)K + ) x B(xco(l-P) -> PP) = (0.030 ± 0.004) x 10" 6 [l5]. We combine our 



upper limit for X(3872) with the known value for B{B + -> X(3872)K + ) x £(X(3872) ->■ 



10 



Table 3: Signal yields, efficiency ratios, ratios of branching fractions and corresponding upper 
limits. 



B + — > (mode) 


Yield 


Cmode/ ^J/V> 


7£(mode) 


Upper Limit 


— > ppK + 


± stat ± syst 


± syst 


± stat ± syst 


95% CL 


J/ip 


1 A f d i i o a 

1458 ± 42 ± 24 




1 




total 


6951 ±176 ±171 


0.970 ±0.002 


4.91 ± 0.19 ± 0.14 




M pP < 2.85 GeV/c 2 


3238 ±122 ±121 


1.097 ±0.006 


2.02 ± 0.10 ± 0.08 




T]c{lS)K+ 


856 ± 46 ± 19 


1.016 ±0.034 


0.578 ±0.035 ±0.026 




i)(2S)K+ 


107 ± 16± 13 


0.921 ±0.044 


0.080 ±0.012 ±0.009 




T) C (2S)K+ 


39 ± 15 ± 5 


0.927 ±0.041 


0.029 ±0.011 ±0.004 


< 0.048 


Xco(lP)K+ 


15 ± 13 ± 4 


0.957 ±0.024 


0.011 ±0.009 ±0.003 


< 0.028 


h c {\P)K+ 


21 ± 11 ± 5 


0.943 ±0.032 


0.015 ±0.008 ±0.004 


< 0.029 


X(3872)K + 


-9± 8± 2 


0.896 ±0.058 


-0.007 ±0.006 ±0.002 


< 0.008 


X(3915)K+ 


13± 17± 5 


0.890 ±0.062 


0.010 ±0.013 ±0.002 


< 0.032 



Table 4: Branching fractions for B + — > (mode) — > ppK + derived using the world average value 
of the B(B + —7- J/ipK + ) and B{Jfy — > pp) branching fractions 15 . The first uncertainties are 



statistical, the second systematic in the present measurement, and the third systematic from the 
uncertainty on the J/tp branching fraction. 



B+ 


B(B + -)■ (mode) ->■ ppK + ) 


UL (95% CL) 


Previous measurements 


decay mode 


(xlO 6 ) 


(xlO 6 ) 


(xlO 6 ) [lo] 


total 


10.81 ± 0.42 ± 0.30 ± 0.49 




10.76^;| ± 0.70 


M pP < 2.85GeV/c 2 


4.46 ± 0.21 ± 0.18 ± 0.20 




5.12 ±0.31 


Vc(lS)K+ 


1.27 ± 0.08 ± 0.05 ± 0.06 




1.54 ±0.16 


i){2S)K + 


0.175 ±0.027 ±0.020 ±0.008 




0.176 ±0.012 


rj c (2S)K+ 


0.063 ±0.025 ±0.009 ±0.003 


< 0.106 






0.024 ± 0.021 ± 0.006 ± 0.001 


< 0.062 


0.030 ±0.004 


h c (lP)K + 


0.034 ±0.018 ±0.008 ±0.002 


< 0.064 




X(3872)K + 


-0.015 ± 0.013 ± 0.003 ± 0.001 


< 0.017 




X(3915)K+ 


0.022 ± 0.029 ± 0.004 ± 0.001 


< 0.071 





J/ipTr+Ti-) = (8.6 ± 0.8) x 10" 6 (15) to obtain the limit 

£(X(3872) pp) 



i3(X(3872) J/iJ7i+n-) 



< 2.0 x 10" 



This limit challenges some of the predictions for the molecular interpretations of the X(3872) 
state and is approaching the range of predictions for a conventional xa(2P) state [l6|[T7] . 
Using our result and the r] c (2S) branching fraction B(B + — > r] c (2S)K + ) x B(i] c (2S) — > 
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KKtc) = (3A±l%) x 10 



15 , a limit of 

B( Vc (2S) -» pp) 
B(rj c (2S) ->• KKtx) 



< 3.1 x 10 



-2 



is obtained. 



8 Summary 

Based on a sample of 6951 ±176 B + — )■ ppK + decays reconstructed in a data sample, 
corresponding to an integrated luminosity of l.Ofb -1 , collected with the LHCb detector, 
the following relative branching fractions are measured 



B(B + ppK + ) 



total 



B{B+ J/i/jK+ ppK+) 

B(B + -> ppK + ) M P p<2.85GeV/ C 2 


B(B- 
B(B + 




— )■ ppK + ) 
v -> ppK + ) 


B[B~ 
B(B + 


-»■ ip(2S)K' 


— >• ppK + ) 
h -> ppK + ) 



B(B+ J/^K+ ppiT- 



4.91 ±0.19 (stat) ±0.14 (syst), 
: 2.02 ± 0.10 (stat) ± 0.08 (syst), 
0.578 ± 0.035 (stat) ± 0.025 (syst), 
0.080 ± 0.012 (stat) ± 0.009 (syst). 



An upper limit on the ratio B ^^^j^k+^pk+) " < 0-017 is obtained, from which a limit 

B(*(3872)->pp) <2Qxl0 -3 
£(X(3872) J/^tt+tt-) 

is derived. 
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